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Abstract 

We present an evaluation of the fluctuation determinant which appears as 
a prefactor in the instanton transition rate for the two-dimensional Abelian 
Higgs model. The corrections are found to change the rate at most by a 
factor of 2 for 0.4 < Mh/Mw < 2.0. 
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1 Introduction 



The Abelian Higgs model in (1 + 1) has found considerable attention recently 
since on the one hand it shares certain features with the electroweak theory 
and on the other hand it is simple enough to serve as a theoretical and 
numerical laboratory. 

In the context of the baryon number violation the high temperature 
sphaleron transition has been studied [|l|-0, for which exact classical solu- 
tions and an exact expression of the sphaleron determinant 0-0] are known, 
providing thus a complete one-loop semiclassical transition rate which can 
be studied numerically on the lattice, e.g. by measuring the fluctuations of 
the Chern-Simons number. 

Another prominent feature of the model is the existence of instanton 
solutions that give rise again to fluctuations in the topological charge of 
the vacuum and thereby to baryon number violation. It has also been used P| 
to study the possibility of baryon number violation in high energy scattering 
processes. 

If the parameters of the model are chosen appropriately, the instantons 
are suffifiently rare and it should be sufficient to consider a dilute gas of 
instantons with Chern-Simons charge q = ±1. For such a dilute gas of 
instantons the transition rate, or equivalently the density of instantons in 
the Euclidean plane, is given by 

r = ^^v-^'^eM-s{<Pci) - ^ct(0d)) (1.1) 

to 1-loop accuracy. Here S{(l)ci) is the instanton action. The coefficient V 
represents the effect of quantum fluctuations around the instanton conflgura- 
tion and arises from the Gaussian approximation to the functional integral. 
This is the object whose computation we will consider here. It is given in 
general form by 

det (M) , i_ioopN r,\ 

^ = d<M^ = ^"P^'^^/^ ^ 

where the second equation relates it to the one-loop effective action. The 
operators are the fluctuation operators obtained by taking the second 
functional derivative of the action at the instanton and vacuum background 
fleld conflgurations. The prime in the determinant implies omitting of the 
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two translation zero modes. The first prefactor S{(f)ci)/'2TT takes into account 
of the integration of the translation mode collective coordinates. Finally the 
counterterm action in the exponent will absorb the ultraviolet divergences 
of V. One may also include a corresponding determinant for fermions, which 
for massless fermions is even known analytically (see below). However in 
lattice simulations the instanton rate and therefore fermion number violation 
can be measured by studying fluctuations of the Chern-Simons number and 
it is not necessary to include fermions. 

For Mh/Mw 7^ 1 even the classical instanton profiles are known only 
numerically, so the evaluation of the effective action has to be performed 
numerically. A method for such computations has been proposed previously 
[ PII ; it has been used recently for the computation of the fluctuation deter- 
minant of the electroweak sphaleron |Tl|] and for the case considered here in 
[ P^ , on which the present work is based. 

This paper is organized as follows: In the next section we outline the basic 
relations of the Abelian Higgs model. The fluctuation operator is derived in 
section 3, its partial wave reduction in section 4. The method of computation 
is presented in section 5. In section 6 we consider the renormalization of the 
effective action and the removal of zero modes. The results are presented 
and discussed in section 7. 



2 Basic relations 

The Abelian Higgs model in (1+1) dimensions is defined by the Lagrange 
density (written in the Euclidean form relevant here) 

^ = \{F,.f + \\D,<P\^ + ^ (l^r - + . (2.1) 

Here 

Uf rif 

i=i j=i 



2 



The particle spectrum consists of Higgs bosons of mass = 2Af^, 
vector bosons of mass = g'^v'^ and left (right) handed massless fermions 
of charge g. The anomaly of the gauge invariant fermionic current 

Uf rif 

j.=i:^fi,^f+i:^Ri.'^'k^ (2-2) 
i=i j=i 

is given by 

d,J, = 2nj {£e,.F,,^ . (2.3) 

The integral over the divergence of the current which measures the baryon 
number violation is given by 

AF = 1 d^xd^J^ (2.4) 

= 2nf (J^ J £xe^,F^,^ = 2nfq , (2.5) 

Here q denotes the Chern-Simons term in two dimensions (siehe z.B. |T3|1 ) 
and baryon number violation is therefore related to euclidean gauge field con- 
figurations with nonvanishing topological charge q. These are the instanton 
solutions which mediate tunneling transitions changing the topological charge 
by q units. We will assume here that the instantons transitions are described 
sufficiently well by a dilute gas of instantons with Chern Simons number 
q = ±1, a situation for which the rate formula given in the introduction is 
supposed to hold. 

A structure that exhibits such a topological charge and satisfies the eu- 
clidean equations of motion is given by the Nielsen-Olesen vortex 0. The 
spherically symmetric ansatz for this solution is given by 

Atix) = '-^A{r) (2.6) 

(f)^^{x) = t;/(r)e*^(^) . (2.7) 

In order to have a purely real Higgs field one performs a gauge transfor- 
mation 

Af, A^-d^if/g 
^L(i?) ^ e^'^-^LiR) (2.8) 
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to obtain the instanton fields in the singular gauge 



At{x) = ^(A(r) + 1) (2.9) 
(f)'\x) = vfir) (2.10) 
With this ansatz the euclidean action takes the form 

H- ^(f('-)-l)^) (2.11) 

For the case = My/ an exact solution to the variational equation is 




known [|T4| for which the classical action takes the value Sd = tcv"^. We will 
consider here the general case, however, for which the classical equations of 
motion 

^^^^-^(fW-l))/(r) = (2.12) 

- — - mj^f (r)l ^(r) = mj,,/' (ri2.13) 

have to be solved numerically. 

Imposing the boundary conditions on the profile functions 

A(r) const ■ , A{r) —1 

/(r) const ■ r , /(r) 1 




(2.14) 



the Chern-Simons number is 1 and the action is finite. 

Since we will consider fluctuations around these solutions later on a good 
numerical accuracy for the profile functions /(r) and A{r) is required. We 
have found that the method used previously by Bais and Primack |jl5[ in 
order to obtain precise profiles for the 't Hooft-Polyakov monopole is very 
suitable also in this context. The values for the classical action - which 
determine also the translation mode prefactor - are given in Table 1. They 
agree almost perfectly with with the results of Jacobs and Rebbi |jl6[. 

The classical action in units of vry^ is plotted in Fig. 1 for .4 < Mh/Mw < 
2. The plot suggests a power law behaviour (Mh/Mw)'^ where 0.40 < p < 
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.43; our data and those of Ref. ITBi are precise enough to rule out an exact 



power dependence. While the best fit is obtained with p ~ .41 a suggestive 
number in this range is p = 1 — 7, where 7 is Eulers constant. This could be 
an asymptotic dependence for large M^/M^y, it is also displayed in Fig.l. 

Though we are not interested here in the effect of fermionic fluctuations, 
we could not resist to use our profiles to calculate the fermion determinant for 
massless fermions. It is known exactly [17 ; expressed as an effective action 
it is, per (left plus right handed) fermion and with fermionic zero modes 
removed, 

gferm = f ^xad'^a (2.15) 

where for the instanton 

a(r) = r dr^ (2.16) 



r 

The results are given in Table 1 and plotted with the bosonic effective action 
in Fig. 4. 



3 The fluctuation operator 

The fluctuation operator is defined in general form as 

where ipi denote the fluctuating fields and ipf the "classical" background field 
configuration which here will be the instanton and the vacuum configurations. 
If the fields are expanded around the background field as ipi = ipf + (pi and 
the Langrange density is expanded accordingly then the fluctuation operator 
is related to the second order Lagrange density via 

= \<P:M.A (3.2) 

In terms of the fluctuation operators M. on the instanton and M.^ on the 
vacuum backgrounds the effective action is defined as 

S,„ = UJ'^\ (3.3) 

2 IdetA^of ^ ' 
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For our specific model we expand as 

A, = A^l^+a, (3.4) 
= <p'' + V, (3.5) 

In order to eliminate the gauge degrees of freedom we introduce, as in 
Ref. 0, the background gauge function 

:F{A) = d,A, + I ((0^') V - <p^'<p*) (3.6) 

which leads in the Feynman background gauge to the gauge fixing Lagrange 
density 

= l{d,a,Y - '-^a,{^d,<P'^' + <P^'d,^ - ^*d,<P^' - <P^'d,^*) (3.7) 

The corresponding Fadeev-Popov Lagrangian becomes 

Cpp = lv*i-d' + g'{ry)v (3.8) 

In terms of the real components = (pi + i(p2 and r] = {r]i + ir]2)/V^ the 
second order Lagrange density becomes now 



2 

+V,l{-d' + g'Al + X{3<l>'-v'))^, 



'2 



+ (p2{gA^dy)ipi + ipi{-gA^dy)ip2 (3.9) 
+ a^,{2g'^A^(j))ipi+ a^,{2gd^(l))ip2 
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Specifying now the fluctuating flelds ((61, 02, 4>3, 04, 05) as (oi, 02, (fii, (p2, V12) 
the nonvanishing components of Ai are 

Mu = -d' + g'<P' 
Mi3 = 29^A,<p 
M23 = 2g^A2(j) 

Mu = -9^ + /4 + V + A(0^ - 

, 2 



A^22 






= 2gd,(j) 


M24 


= 25(920 


M34 




M43 





It is understood that the contribution of the Fadeev-Popov operator M55 
enters with a negative sign and a factor 2 into the definition of the effective 
action. The fluctuation operators for the instanton and vacuum background 
are now obtained by substituting the corresponding classical flclds. The 
vacuum fluctuation operator becomes a diagonal matrix of Klein-Gordon op- 
erators with masses {Mw, Mw, Mw, Mh, Mw)- It is convenient to introduce 
a potential V via 

M^M'^ + V (3.10) 
The potential V will be specified below after partial wave decomposition. 

4 Partial wave decomposition 

The fluctuation operator Ai can be decomposed into partial waves and its 
the determinant decomposes accordingly. 

+00 

lndetA<= E IndetM^ (4.1) 

n=— 00 

We introduce the following partial wave decomposition for fields 

+°° ^ cos A e^"'^ . . - sm(p\ e'""^ 



+00 



-00 pin(p 
n=-oo x/27r 



-00 



<P2 ^ ^n{r) 



27r 
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+0O pintp 

Vl2 = Yl 9n{r) 



2ti 

After inserting these expressions into the Lagrange density and using the 
reahty conditions for the fields one finds that the following combinations are 
relatively real and make the fluctuation operators symmetric: 



Ff(r) 




Cn(r)) 






Cn(r)) 




= K{r) 






= iK{r) 






= 9n{r) 





Writing the partial fluctuation operators - omitting the index n in the fol- 
lowing - as 

M = M° + V (4.2) 
the free operators M*^ become diagonal matrices with elements 

M° = --^-i-^ + 4 + M^ (4.3) 

where (rij) = (n — + n, n) and (Mj) = {Mw, Mw, Mw, Mh, Mw)- 
The potentials V takes the elements 



r- A+1 

r 

^ 13 

2 ^2 





= rnl,{p- 


-1) 


12 




= V2m^f 




14 


* 22 


= vi^, 




^ 23 








^33 






n 




''55 




-1) 


V,5 



= 



Chosing the dimensionless variable Mwr one realizes that the fluctuation op- 
erator depends only on the ratio Mh/Mw up to an overall factor which 
cancels in the ratio with the free operator. 
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5 Computation of the fluctuation determi- 
nant 



The method for computing the fluctuation determinant used here is based 
on the use uf the Euchdean Green function of the fluctuation operator. This 
Green function is defined by 

{M + u^)g{x,x\u)^l5{x-x'). (5.1) 

and similarly for the operator Ai^. It contains the information on the eigen- 
values Aq, of the fluctuation operator via 



f£xTvg{x,x,u)^J2-^. (5.2) 

If we define the function F{y) via 

F{v) = j (fxTr{g{x, X, u) - g°{x, x, u)) (5.3) 

we have 

-I "-'^^•'^^^Hii^T^fK^)} 

For e — > this is just the logarithm of the ratio of fluctuation determinants, i. 
e. the one loop effective action, regularized with a Pauli-Villars cutoff. The 
regularization can be removed, the integral can be taken to infinity, after 
subtracting the one loop counterterm action (sec below). Before taking the 
limit e — > the two zero eigenvalues have to be removed by subtracting their 
contribution Ine^. Of course e has the dimension of energy. We used here 
the scale Mw throughout, i.e. by making the radial variable dimensionless. 
So Seff contains now a term — InM^, the numerical prefactor and 
therefore the rate are computed in units of M^. 

After these more formal considerations we have to present a practical way 
for computing F(i/). This is done by using the partial wave decomposition 
to write 

+ 00 

F{iy)= Y: Fn{y) (5.5) 

n=—oo 

where 

F„(z/) = j drrTr(G„(r,r,z/) - G°(r,r,z/)), (5.6) 



9 



the partial wave Green functions being defined by 

(M„ + iy^)Gn{r, r', u) = l-6{r - r') (5.7) 

r 

For the Green function is simply by a diagonal matrix with elements 

G° ,,ir,r',u) = I^X^,r^)K^XK,r^) (5.8) 

where k = \J Mf + z/^. For the Green function of the operator M„ the matrix 
elements become similarly 

G..,(r,r',i.) = /r(r<)/;+(r>) (5.9) 

where the functions /f"^ form a fundamental system of solutions of (|5.71 ) 
regular as r — >■ for the minus sign and as r — > cxd for the plus sign. The 
correct normalization is obtained by imposing the boundary conditions 

ftHr) ^ 6:K^X^,r) (5.10) 

as r oo. Actually we have solved numerically the differential equations for 
the functions h'^^ defined by 

ft^ = B^X^,r){6r + ht{r)) (5.11) 

where are the appropriate Bessel functions. In this way one keeps track 
of the free contribution oc 6f and one has 

Tr (G(r,r,z/) - G°(r,r,z/)) 

= (hlir) + hl^ir) + hr{r)hr{r))I^Mr)KnX^^r) (5.12) 

to be inserted into (|5.5|) . The partial wave contributions behave as for 
large n. The summation implied by Eq. ( |5.5| ) has been performed up to 
maximally n = 25, the asymptotic tail was appended by fitting the last five 
terms with an expression = c^n'^ + C4n~^ + c^n'^ and adding the sum 
over the a„ from n + 1 to oo. The convergence was monitored by applying 
this procedure already in each step of the n summation taking n = n. The 
convergence was found to be excellent up to values of v of the order 5. It 
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has to be said, though, that there is considerable cancellation between the 
negative n = contribution and the higher terms. Indeed the v integration 
over the n = term alone would be divergent even after renormalization. 
This seems to be an inherent feature for functional determinants for topolog- 
ically nontrivial configurations, it is related to the fact that the centrifugal 
barriers of the operator M„ at r = are different from those of M°. This 
feature makes also a direct application of a theorem on functional determi- 
nants used for the faster and more elegant method of Ref. [jl8[ impossible. 



The deformation of the centrifugal barriers is not related to our using the 
singular gauge for he classical instanton field. In fact it can be shown by di- 
rect calculation that the fluctation equations do not change under the gauge 
gauge transformation (^.81). 

Fortunately the asymptotic behaviour oi F{v) which is as after renor- 
malization sets already in when this function has dropped to values of order 
10~^ and there the cancellation is not yet delicate. 

There is a further problem which we have to address here which is related 
to the coupling of fields with different mass in the system of gauge, Higgs 
and Goldstone fields. While normally the solutions of the coupled system 
fulfil vacuum boundary conditions at r ^ oo, i. e. the potential decreases 
sufficiently fast, the cross terms Vj4 can cause the Higgs field to change 
the asymptotic behaviour of the gauge and Goldstone components. The 
solution regular at r = behaves normally as exp(Kjr). If the physical Higgs 
component is multiplied by Vj4 one obtains a behaviour exp((Kiy — Mw)^). 
This expression enters the right hand sides of the equations for the Goldstone 
and gauge fields which themselves behave as exp(Ki4/T'). So if Mh > IMw 
these fields change their asymptotic behaviour. We find that the integral of 
the trace of the Green function over r ceases to exist. We think that this not 
a shortcome of the method but a systematical property of the fluctuation 
determinant. Indeed for Mh > 2Mw the Higgs boson can decay into pairs 
of gauge particles and also the singularity structure of perturbative graphs 
changes qualitatively. This subjects merits further consideration; here we just 
restrict our computation to Higgs masses smaller than 2Mw- The gauge fields 
cannot, on the other hand, decay into Higgs particles, since their coupling 
joins a physical and a Goldstone Higgs; indeed our coupled system has no 
problems of principle for small Higgs masses. 
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6 Renormalization and zero modes 



The Abelian Higgs model is super-renormalizable; all divergences can be 
removed by a mass counterterm for the Higgs field and a counter term for 
the vacuum loops. Expanding around (j) = v and using the corresponding 
Feynman rules we find divergent tadpole diagrams of the form represented 
in Fig. 2, where the internal lines are the various Higgs, vector and Fadeev- 
Popov fields. The various couplings can be read off from the second order 
Lagrangean ( |3.9| ). For the vertices of the second graph we find — 3iA/2 for the 
physical Higgs of mass M^, —i{g'^ + A)/2 for the Goldstone Higgs of mass 
Mw, ig^g^J.u/'^ for the gauge field and —ig^/2 for the Fadeev-Popov fields. 
For the first graph we find the same vertex factors multiplied by 2v. As 
a consequence, in summing the contributions from both graphs the external 
line factors combine as {(p—vY +2v{(j)—v) = {(p'^—v'^). The contributions from 
the gauge field and Fadeev-Popov loops cancel as they should. The tadpole 
graphs with external gauge field lines (not presented in Fig. 1) cancel against 
second order graphs as usual in scalar QED. The counter term action takes 
the form 

S,t = Um' j £x{<P^-v^) (6.1) 
where in unregularized form 

<fk 1 , . r d^k 



5m' = 3A 



(27r)2A;2 + M^ 'J (27r)2 P + M^, 

In the Pauli-Villars regularization chosen her we rewrite the divergent mo- 
mentum integrals via 

d^k f 1 1 \ 



(27r)2 \k^ + Mf k^ + Mf + Ay 



(27r)2 Jo (p + M2 + z/2)2 

1 1 



271 Jo M2 + Z/2- ^^-^^ 

so that the divergent terms can be rewritten directly as a contribution to a 
counterterm -Fcj(z/) in the integral over u. We find 
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If Fctiiy) is subtracted, Fren{i^) = F{v) — Fct{i^) behaves as u^^ as u ^ oo 
and the Pauh-Villars cutoff, i.e. the upper limit of integration, can be sent 
to oo. 

Instead of subtracting the tadpole contributions from F{i>) this subtrac- 
tion can be performed already in the partial waves. The tadpole terms can 
be easily recognized in the potential given at the end of section 4 as the 
diagonal terms proportional to (/^ ~ !)• Denoting these terms by V*""^ their 
contribution to the first order Green function becomes 

Gniir: r\ u)^J drVGl ,,(r, r", u{r": r\ v) (6.5) 

(no summation over i). Using some identities for Bessel functions it can be 
shown that after taking the trace, integrating over r and summing up the 
partial waves one obtains Fdiy)- In the actual computation we have removed 
the tadpole contributions directly in the partial waves. As an illustration 
we show however, in Fig. 3, the function F{y) before the subtraction of 
the tadpole and Fadeev-Popov contributions, both of these contributions, 
and the final F'^^'^iy). It follows from perturbation theory that the former 
behave as v""^ asymptotically, while the latter behaves as v~^. The numerical 
integration was performed up to to region where the asymptotic behaviour 
sets in. The remaining integral was performed as / dvv-^ with a coefficient 
determined by the last point. The contribution of the integral from v^nax 
to oo is of the order .05 and the error introduced by the extrapolation is 
certainly one order of magnitude smaller than this value. 

One notes in Fig. 3 that F{v) behaves for small z/ as a behaviour 

that is due to the translation zero modes and makes the subtraction of In 
necessary while the lower limit of the integration is taken to 0. In practice, 
the zero mode pole appears slightly shifted to = Aq ~ .02 as can be seen 
from the departure of the expected behaviour for i/ < .1. So a term ln(e^ — Aq) 
has to be subtracted instead. The integrand was, for u < 1 decomposed into 
a pole term and a finite contribution and the former one was integrated 
analytically. Aq can be fixed to at least three significant digits and the finite 
term turns out to show a smooth behaviour oc i/; we think that this procedure 
introduces an error of Seff below .01. So including the estimate for the error 
in the asymptotic extrapolation and another .05 (i.e. ~ 10%) for errors in 
the numerical integration we think that we have determined Seff to within 
an error of .07. 
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7 Discussion and Conclusion 



The results of our computation of the one-loop effective action 

Seff = -lim( / dvvFren{^) + Ine^) (7.1) 

arc shown in Fig. 4. The fluctuation prcfactor (including the coun- 

terterm action) is given by exp(— S'^jj). Due to the subtraction of the zero 
mode contribution Ine^ it has dimension {length)^"^ . Since we have used in 
our computation the dimensionless variable Mwt the units for the rate are 
(the action and therefore the zero mode prefactor being dimensionless) . 
As mentioned in the previous section we estimate the error of our numerical 
result for Seff to be of the order of .07 units. 

In contrast to an analoguos computation of the fluctuation prefactor for 
the sphaleron transition in the electroweak theory here the effects of the 
quantum fluctuations on the transition rate remain quite small, less than a 
factor of 2. This could have been expected on the grounds that the number 
of fluctuating fields small; effectively - in view of the cancellation between 
gauge field and Fadccv-Popov degrees of freedom - we have just the physical 
and the Goldstone part of the Higgs field. Also the dimension of space is 
reduced from three to two. We think nevertheless that this expectation had 
to be checked by a direct computation. 

One cannot compare the classical and quantum action without specifying 
the dimensionless vacuum expectation value v = Mw/g"^- If f — 1 then the 
classical action is ~ tt x 0(1). This has to be considered as an absolute 
lower limit if one wants to justify the dilute instanton gas approximation. 
The fact that the one loop correction is then one order of magnitude smaller 
supports the use of the semiclassical approximation. It would be interesting 
to compare it to lattice simulations. 
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Figure Captions 



Fig. 1 The classical action. We present the action in (dimensionless) units 
TTf^. The full circles are the numerical results, the curve displays a 
possible asymptotic dependence (M^/Mvy)*^! — 7), where 7 is Eulers 
constant. 

Fig. 2 The divergent tadpole graphs. The vertex factors and internal line 
masses are given in the text. We do not display the analoguous graphs 
with external classical gauge field legs since they cancel against second 
order contributions. 

Fig. 3 The integrand -F(i^) for Mh = M^. Empty diamonds: F{i>) before 
tadpole and Fadeev-Popov subtraction; empty circles: Fct{i'); empty 
squares: the Fadeev Popov term; full circles: Fren{j^)- The dotted hne 
corresponds to the behaviour due to the zero modes. The dashed 
line shows the extrapolated asymptotic z/~^ behaviour. 

Fig. 4 The one loop effective action. The vertical lines are the numerical 
results for the bosonic effective action Sljj-""^, the length of the lines 

indicate the error; the empty squares are the effective action Sl^P for 
massless fermions given by Eq. (|2.15|) . The curves are spline fits. 
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Tables 





Sh 




q1 — I()(>P 


0.40 


0.696196 


-0.37853 


0.315 


0.60 


0.813053 


-0.38459 


0.233 


0.80 


0.912305 


-0.39051 


0.102 


1.00 


1.000000 


-0.39616 


-0.025 


1.25 


1.097914 


-0.40277 


-0.156 


1.50 


1.186013 


-0.40886 


-0.276 


1.75 


1.266416 


-0.41445 


-0.379 


2.00 


1.340550 


-0.41956 


-0.461 



Table 1: Classical and one-loop actions for various values of Mh/Mw- (left 
plus right-handed) fermion. S],Jf°°^ is the one loop bosonic action computed 
here. 
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